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ABSTRACT 

We investigate the decay mechanism in the — > (f)(j)K^ decay with the 00 invariant 
mass below the charm threshold and in the neighborhood of the r]c invariant mass region. 
Our approach is based on the use of factorization model and the knowledge of matrix 
elements of the weak currents. For the B meson weak transition we apply form factor for- 
malism, while for the light mesons we use effective weak and strong Lagrangians. We find 
that the dominant contributions to the branching ratio come from the 77, r]' and ?7(1490)) 
pole terms of the penguin operators in the decay chains B~ —>■ ri{r]', ?7(1490))i^'~ — *• (jxpK' . 
Our prediction for the branching ratio is in agreement with the Belle's result. 



I. INTRODUCTION 



It is a very fruitful era in B meson physics. A lot of experimental data on B meson 
decays is coming from the B meson factories. Many of their results are still not explained. 
Recently, Belle collaboration has announced the observation of the BR(i?^ (fxpK^) = 
{2.6tli ± 0-3) X 10"^ [T for a 00 invariant mass below 2.85 GeV. This is the first of the 
three-body B decays with two vector mesons and one pseudoscalar meson in the final 
state that has been observed. The B meson decays into three pseudoscalar mesons have 
been studied [21 El within heavy quark symmetry accompanied by chiral symmetry. One 
might explain the observed rates using heavy quark symmetry for the strong vertices, 
while for the weak transition we rely on the existing knowledge of the form factors 2 . 
The three-body decay with two vector meson states and one pseudoscalar is much more 
difficult to approach. 

The additional insight on the decay mechanism might come from the analysis of the 
B meson two-body decays. Particularly interesting are the decays B^ B^ 
ri{r]')K^ and B K*^(j). They have been extensively studied using different existing 
techniques: the naive factorization [U El E], the QCD factorization [71 and the SU(3) 
symmetry Each of these decay modes is rather difficult to explain theoretically. The 
decays B^ — >■ (j)K^ and B^ — ^ K*^(f) might have significant annihilation contribution 
lUElE], but it is not simple to have a consistent treatment of it. There is an interesting 
proposal [HI in which the angular distributions of the final outgoing particles can be used to 
estimate the magnitude of annihilation contribution to the amplitude. However, we have 
to wait for the new experimental data to extract the size of the annihilation contribution. 
The B^ —y r]{ri')K^ decay rate has not been easy to explain. It accounts for the well- 
known problem of the t] — rj' mixing piU ITT] as can be seen from a variety of approaches 
used for B^ r]{r]')K^ jUH^lilSllIll- In the B^ r]{r]')K^ decay mode, it seems that 
the annihilation contribution is not very significant |H |^ . 

One has to expect that the above described difficulties in these decay modes might 
appear in the three-body decay we discuss. Based on the current knowledge of two- 
body transitions, we build a simple model which might clarify the role of the non-charm 
contributions in the BR(i?^ (fxpK^) decay. In our study of the B^ —>■ cfxpK^ decay 
mechanism, we follow the assumption in ref. |2| and use double and single pole form 
factors for the B meson semileptonic transitions [T^ lltij. Our approach is based on naive 
factorization, as QCD factorization has not been developed yet for three-body decays. The 
SU{?)) symmetry approach is not applicable due to the limited number of the observed 
B decay modes. In our model we keep only dominant contributions and as in the case 
of two-body charmless B decays, we do not include annihilation contributions. We use a 
pole model including the low-lying meson resonances and possible contributions coming 
from higher mass excited states. 

In order to compare our result with the Belle's result, we include in our calculation 
the interference between the non-resonant 5" K^(j)(j) and the resonant B^ K^rjc —>■ 
K~(l)(f) decay amplitude. 

In Section II we present the basic elements of our model, while in Section III we give 
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the results for the three-body decay amplitude and discuss possible contributions to the 
decay rate. 



II. THE MODEL 

The ub suss transition which can produce two (f) mesons in the final state via strong 
interactions, can be realized by the effective weak Hamiltonian [T7]-|2Uj: 
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where Oi and O2 are the tree-level operators, O3 — Oq are gluonic penguin operators 
and O7 — Oio are electroweak penguin operators. Superscripts u, c, t on the Wilson co- 
efficients denote the internal quark in penguin loop. In order to apply the factorization 
approximation we rearrange the above operators using Fierz transformations and leave 
only color-singlet ones. One then comes to the effective weak Hamiltonian given by Eq. 
([Q) replacing the coefficients q by a^. The relevant operators are: 



{ub)v~A{su)v-A, 
{sb)v-A{uu)v-A, 
E,Ol = E,isb)v-A{qq) 
T.,Ol = Y.,{sq)v-Am 
EaOl = Eai^b)v.Aiqq) 



V-A, 
V-A, 
V+A, 



= -2 01 = -2 + 75)g)(g(i - 75)^) 

Or 

0,0 



E,0'7 = E,le,{sb)v-Aiqqh 



— 2Eje,(.-(i+75)g)(g(i-75)&), 

Eg C'g = Eg leq{sb)v-A{qq)v-A, 
q 2*^ 



IV-A 



E,0% = E,leq{sq) 
The Wilson coefficients are taken from [TUj : 



V-A, 



ai = 1.05, 

04 = (-43 - 16i) X 10-3, 
07 = (0.4-0.9z) X 10-^ 



02 = 0.07, 
ag = -53 X 10-^ 
as = (3.3 - 0.3i) x 10~\ 
= (-13-0.3i) X 10-^ 



as = 47 X 10-^ 

ag = (-54 - I6i) X 10- 

ag = (-91 - 0.9z) X 10" 



(3) 



For the CKM matrix elements (Vij) we use Wolfenstein parametrization:V^6l^* = — AA^ 
and VubV*^ = v4A^(p - if]), where A = 0.83, A = 0.222, p = 0.217/(1 - AV2) and 
fj = 0.331/(1 — A^/2). The standard decomposition of the weak current matrix elements 
is: ^ 

„7 \q ) , D- ^ ' u, A t 2\ 



{V{k,e,my)\qV^q\P{p,M)) 



_2VV_ 

M + my 



+ 2imv^q^Ao{ 



+i{M + mv] 



g2 



Mq' 



e-q 



M + mv 



M2 



2 



A2{q' 



(4) 



form factor 


F+ 




V 


^0 




A2 


/(O) 


0.36 


0.36 


0.44 


0.45 


0.36 


0.32 


0"! 


0.43 


0.70 


0.45 


0.46 


0.64 


1.23 




0.0 


0.27 


0.0 


0.0 


0.36 


0.38 


M [GeV] 


5.42 


5.42 


5.42 


5.37 


5.42 


5.42 



Table 1: The B ^ K, K* form factors at = and the pole parameters jl6j . 



{Pik,mp)\qT>^q\Pip,M)) = - F^iq') + ^-^^^^^q^Foiq') , (5) 

where q^^ = — and = p^^ + k'^ . Also 

(P(p)|g7'^(l - 75)g|0) = tfpp^ , (0|g7^g|l^(p)) = (?y£^ . (6) 

Using experimental data the decay constants are found to be {g^pl = 0.24 GeV^, 
Igxl = 0.19 GeV^, /k = 0.16 GeV and U = 0.132 GeV. The lattice calculation |22] gives 
for the B meson decay constants J'b = 0.173 GeV and fs^ = 1.22/^. We also take 
gs* = Mb'/b IS]- The dependence of the form factors is studied in [12], where a 
quark model is combined with a fit to lattice and experimental data. This approach 
results in a double pole dependence of F+{q'^), V{q'^) and ^0(9^) 

f(q^) = m (7) 

while for 2(5^) and -Fo('?^) [IH] we have: 

m 



fid' 



[1 - (TgVM2 + a^q^/M^) 



Values of M, /(O), ai and (T2 are listed in Table 1. 

In the evaluation of Oq operator we have as usual ^3] 

— ^ 

^17552 = ■ d^{qiYl5q2) , (9) 

nil + "^2 

9ig2= ^^(917^92). (10) 

mi — m2 

The effects of strong interactions of light mesons are taken into account by using the 
following effective Lagrangian (2211211123: 

Strong = ^Tr(p^[n, d,U]) - i^XXILef^^-f^Trid.p^d^ppU) , (11) 



where 11 and are 3x3 matrices containing pseudoscalar and vector meson field operators 
respectively and / is a pseudoscalar meson decay constant as in Eq. We take 

Cvvn = 0.31 In order to include SU(3) fiavor symmetry breaking, instead of the 

coupling constant coming from the p nn decay {gpmr = 5.9), we use the coupling 
constant coming from the (p — > KK decay rate. Thus, we have g^xK = 6.4. 

For the description of strong interactions between heavy and light mesons, we use 
definitions given in ^H] and heavy quark effective theory to get: 

< Kip,)B:ip,,6)\B{p^+p2) >= ^^{pi + 2p,)^e^ , (12) 
< 4>{Pi,£i)B*{p2,e2)\B;{pi +P2,e) >= ^gBsBs<t,{Pi ■ 6261 ■ e - pi ■ e 62 ■ si) , (13) 

where 5'B.B.</./B./2m^ = 1.5 ± 0.1 and gBiBKfBs/'^^Bt = 0.65 ± 0.05 HH]. 

To account for the rj — rj' mixing, we follow the approach in . Using the quark basis 
(?7g ~ [uu + dd) I \/2 and r]s — ss) , the mixing is given by 

with the mixing angle = 39.3° ± 1.0°. The 77, 77' decay constants are defined by 

< r]\qY{i - 75)g|o >= tf^ , < v'lqYi^ - 75)g|0 >= tf^, , (15) 

where 

= fu,d COS (j)/V2 , = sin , 
/^'''^ = fu,d sin (f)/V2 , = fs cos , 

with fu,d = (1.07 ± 0.02)/^ and /, = (1.34 ± 0.06)/^. The form factors for the B 
transition can be written as: 

^oV(9') =i^o^(?')cos(/)/y2, FS;4q')=F^{q') sin (17) 

The dependence of Fq is described by Eq. (jSI), with Fq{0) = 0.29, (Xi = 0.76 and 
(72 = 0.28 while the dependence of F^ is described by Eq. ((7j), with F^{0) = 0.29 and 
ai = 0.48 HE]. 

Before we consider the B~ —>■ (j)(j)K~ decay amplitude, we check how our model works 
for the two-body decays: B~ —>■ r]K~ , B~ ^'K~ , B~ — > (f)K~ , and B~ (j)K*~ . 
Namely, B~ — > (j)(f)K~ can occur through one of these decay chains: B^ r]{ri')K~ 
followed by 77(77') 00; B^ 4>K~ , followed by K~ K~(f) and B^ (pK*^ followed 
by K*^ — > K^cf). These decays have already been studied within the factorization ap- 
proximation by Ali et. al. Using their formulas for the amplitudes with the Wilson 
coefficients, the form factors and other parameters as given above, we obtain the branch- 
ing ratios for the two-body decays presented in Table 2 together with the experimental 
results. We point out that references jHE], as well as our predictions, include the axial 
anomaly contribution in 6 — > sgg sri{ri'). In our calculation, the contribution of the cc 
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Belle [27. 28j 


BaBar [29, 30, 31, 32 


our model 


B- r]K- 


(5.3li:^ ±0.6) X 10-** 


(2.8l[]:? ± 0.2) X 10"^ 


2.1 X 10-6 


B- r]'K- 




(7.67 ±0.35 ±0.44) x 10"^ 


3.0 X 10-^ 


B- (pK- 


(9.4 ± 1.1 ±0.7) X 10-6 


(10+^1 ±0.5) X 10-6 


14.9 X 10-6 


B- (pK*- 




(12.7+^:^ ± 1.1) X 10-6 


8.6 X 10-6 



Table 2: The experimental and theoretical results for the relevant B two-body decay 
rates. The rates for B~ riK~ and B~ — > ri'K~ are calculated without cc contributions. 
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Figure 1: Feynman diagrams for B (pcpK . 



component in 77,77' was found to be small and therefore we safely neglect it. 
III. THE ^ (p(pK- DECAYS 

The dominant contributions in the B^{p) -^^^(^1)0(^2)0(^3) decay amplitude with 
the 00 invariant mass in the region below the charm threshold are shown in Fig. 1. We 
write the amplitude for this decay in the following form: 

M=aJj2'^i)+^4Y1 ) + -^3^2^ ± A,Ci + {A, + Ae + B)C^ . (18) 

\q=s,u,d / \q=s,u,d / 



Here 



G 

Ciu = (KfeK>2 - VtbVtlias - as ± ag - 07)) 
CI = -^VMa, -a,- l/2{a, - a,))^j- , 
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G ( ( V \ 

Cl^ = ^VtbVtl I I as + a4 - as - l/2{ag + - aj) H —{a& - l/2as) ) -j-- 

V 2 V \ Tribiris J jn 

2 ru 



,a6-l/2a8)^ , (19) 

= ^VubV*^ai - VtbV*{a4 + + (ag + ag)) , 

V2 mbTUs 

Cf = ^(-\4,V;:) (as + a4 + a5 - l/2(a7 + ag + aio)) , (20) 
v2 

where and px are the momenta of t] and meson respectively. The formulas for t]' 
are obtained by replacing and kg with and k'^. The constant fc^ (fc^) projects the ss 
component of the rj [i]') meson and it is equal — sin for rj and cos for rj'. The coefficient 
contains the effect of the axial anomaly as in 01^^. The amplitudes are determined 

by 

Ai = 8iCvvnksFo{x) -e^.^fssi^e'^p'^pl , (21) 

— X 

A2 = SzCvvuhFoix) , _ ^ e^.„^e^e3>^pg , (22) 

fTl^, X 



A, = 8zCyyn^{FS{ml)^^ ^ + (m^))e^.,^4e»f , (23) 



A, = StCvvu-f{F^\ml)^, ^ + (m|))e^.,^e^e3>°pf , (24) 

A o 17 1 2^ Pi •g2P-g3 , p-e2Pi-£3 \ f^^. 

A, = 2V2g,g,.Mm,) _ ^ + , ^ ^ _ _ 2^2 j ^ (25) 

A = -T^?^^^ f 2.(M + m.^rMml) ( w^^^sP?^^ ^ (36) 



— 2m| + — m|^* — x — y 

etiuaf3e2PlP2P§{Pl + P2) ■ £3 , eM'^«/3^3PlP2P3 (Pl + Ps) " ^2 



+\^(m' 



^'k* ~ y + 2^«^ + — ?«|;* — X — y 

(M^ - m|) (m^ - m^) + ^(M^ + 2m^ + m|. - 2x) - 
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<p 



+ 2m^ + m|- — m|-* — x — y 



£2 ■ £3 



M^ + 3ml-x-y ml - m% + y 

X + y — — 2m^ — mj^ + mj^^ mj^* — y 



2pi • £2Pi • £3 



y — A'P — ml X + y — ml — 



+2(M^ + m|i- - a; 



+ 2m^ + — mj^, — x — y m"^* — y 

P3-£2Pl-£3 Pi -£2^2 -£3 



+ 2m| + m|- — m|-« — x — y m\* — y 



2{y -ml + m\)p ■ £2^2 ■ £3 _^ 2(a; + |/ - - - 2mj^)p3 ■ £2 P ■ £3 
mx* — y M"^ + 2ml ~ ~ ^ ~ U 

e - 4M?(~Ki$'-.) ■ ((^^ - "5.)- - Ml(M^ + 4™j - ^D) (27) 

+(M2 - + m^)(p3 • £2P • £3 +P • £2P2 • £3) + {Ml + M^- m^)(p3 • £2^1 • £3 +Pi • £2^2 • £3)) 

In our expressions the two meson polarization vectors are denoted by £2 = £2(^2) 
and S3 = £3(^3), the B~ and B*~ masses are M and Ms, and mjv stands for the mass of 
N meson. 

To obtain the decay width, we make the following integration over the Dalitz plot: 

where y = m\2 = (pi +^2)^, x = m^^ = (p2 +^3)^- Note that we include the factor 1/2 
due to two identical mesons in the final state. In the above integral, upper and lower 
bounds for y are: 



yma.^{E* + E;y-(y/Ef-ml-^/Ef-ml] , (29) 

2 



Vmin = (El + Elf - {^Ef - ml + yjEf - ml^ 



(30) 



with the energies E^ and El given by: 

El^^/2, El^{M^-x-ml)/{2V^). (31) 

The integration over x is bounded by Xmin = 4m^ and Xmax = {M — mK^. 

First we consider only the phase space region with the 00 invariant mass below 
the ric threshold by taking x < (2.85 GeV)^. The Belle collaboration has measured 
BR{B- i^~00)a:<(2.85GeV)2 = {'^■^-qI ± 0.3) X 10"^ while our model gives BR(S- 
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Figure 2: rjcjycf) interaction. 



i^~00)a;<(2,85GeV)2 = 1-8 X 10^^. The Calculated decay rate is the total contributions 
from the parity violating (the terms in amplitude containing e^^ap) and parity conserving 
parts, which do not interfere. The parity violating component gives the rate 1.5 x 10~^, 
while from the parity conserving part we get 0.3 x 10~^. We note that the dominant 
contribution comes from the 77, rj' intermediate states in the graph A1-2 of Fig. 1 and its 
contribution alone gives the branching ratio of 1.3 x 10^®. Since for the B ri{r]')K decay 
rates the annihilation term is not very large we do not expect a significant change in 
the B —>■ (f)(f)K decay rate if its effects are taken into account. 

In addition to the low-lying mesons such as r] and 77', one could expect that higher 
mass excited states in the 1 — 2 GeV region could also make important contribution 
to the amplitude. If the rj.rj' in the diagram, A1-2 (Fig. 1) are replaced by scalar or 
tensor mesons, which contain ss (e.g. /o(980), /2(1270)) one finds that both contributions 
are suppressed. The observed rate 5" /o(980)fC~ [3^ is by an order of magnitude 
smaller than the rate of i? — >■ r]'K~ and the decays of B into a pseudoscalar and a tensor 
meson, are expected to have branching rations of the order 10~^ [21] • The products 
< fo,2\{ub)v±A,s±p\B^ >< K~\{su)v±A,s±p\^ >, (y ± A stands for the left and right 
handed currents, and S ±P are scalar and pseudoscalar densities) can be safely neglected 
because of the small values of the B ^ S,T transition form factor involved in the graphs 
like those in A3-4, (Fig. 1) [SHI- The same arguments hold for the higher mass us excited 
states. 

However, a large contributions to the decay rate can be expected from the higher mass 
excited states with the quantum numbers of 77,77': 7/(1260), 7/(1490) In the study of 
[TT] it has been found that ?7(1295) is most likely (uu + dd)/\/2, while 7/(1490) is almost 
pure ss state. Therefore, we might expect the presence of 7/(1490) in the diagram Ai-2- 
Unfortunately, its interactions are very poorly known and we can make only a very rough 
estimation of the 7/(1490)00 coupling within a naive quark model. The coupling of the rj 
or any state with the same quantum numbers like r/(1490) — 00 could be estimated by a 
ss quark loop triangle graph as shown in Fig. 2. 
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We then find that 



C-n^iq^) (X dx dy{ml-xml-yml + x^ml + y^ml-xy{q'^ -2ml)) K (^2) 
Jo Jo 

Taking the dynamical s quark mass ~ 500 MeV, we roughly estimate Cr^^^im^) : 
Cr,'</,0(m^,) : C'^(i490)<}i</.("^.^(i49o)) = 1 '■ 0-85 : 0.40. We fix Cr,^^{m^2) to be equal to the 
vector-vector-pseudoscalar coupling Cyvp = 0.31. 

Including the contribution of ?7(1490) in the graph A1-2 with above couplings, and 
assuming that there is no axial anomaly term in the coefficient C^^, we find: 

BR{B- ^ ir-00).<(2.85GeV)2 = 3.7 X lO^^ . (33) 

The distribution {l/r)dr/dx as the function of the (jxp invariant mass in the region below 
charm threshold and the Dalitz plot are is given in Fig. 2 for this case. 

Note that the non-resonant contribution in the branching ratio, measured by Belle 
collaboration contains not only the non-resonant amplitude itself, but also the interference 
terms with the resonant contribution K^rjc — ^ K~(j)(f) as in |^. In addition to 

the T]c state there are a number of other cc bound states which might contribute. From 
these, the biggest contribution will arise from the Xco state as its mass is closest to the 
region we discuss {x < 2.85 GeV). This contribution can be obtained from the measured 
B^ — > XcoK~ decay rate One might then expect that the B~ — > XcoK~ (l)(l)K~ 
transition can give additional interference with the calculated rates. However, the rate for 
Xco 00 is ten times smaller than the rate r/c (jxp ^"^^ 'we expect additional suppression. 
This leads us to the conclusion that the interference of the non-resonant and the resonant 
terms from the cc states other than r]c is negligible below the charm threshold. 

Next, we comment on the interference of the r]c resonance with the non-resonant 
contribution in the region of the phase space with the invariant mass of the 00 state 
within the region (2.94 GeV)^ < x < (3.02 GeV)^. The decay rate for B^ K^rjc is not 
theoretically very well understood. Naive factorization leads to a decay rate ten times 
smaller than the branching ratios 6.9^3 q x 10~^ measured by CLEO collaboration j^U], 
(1.34±0.09±0.13±0.41)xl0-^ by BaBar collaboration gT] or (1.25±0.14+°:J°±0.38)xl0-3 
by Belle collaboration QCD factorization seems to face similar problem in explaining 
this decay amplitude jlHl- On the other hand, the decay rjc —>■ 00 rate is not very well 
understood. First, the statistics for the r/c — 00 decay rate is rather poor (the error 
stated in [211 seems to be underestimated jS]). Secondly, by assuming the SU(3) fiavor 
symmetry one cannot reproduce both the rjc — *■ 00 and the rjc —>■ pp measured decay rates. 

Facing these difficulties we use experimental data to estimate the size of this resonant 
contribution. In the phase space region (2.94 GeV)^ < x < (3.02 GeV)^, Belle measures 
BR(5- ^ K-T]^) X BR(r/c ^ 00) = (2.2+^;° ± 0.5) x lO'^. We model BR(B- 
K~r], T]') X BR(?7, rj' 00) by taking the rjc propagator and by fitting the Belle's data. 
The results for the interference are given on Fig. 3, where we present both cases: positive 
and negative interference terms, which is the result of an unknown phase in the rjc 00 
decay amplitude. 
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dr/rdx [10" 

. 08r 




x[Gev^ ] 



17.5 

15 
12 . 5 
^ 10 
7.5 



2.5 




Figure 3: spectrum for the B K (fxj) decay with the 00 invariant mass in the 
region below charm threshold and the Dalitz plot. 
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dr/rdx [10"^] 




Figure 4: The spectrum for B~ K'cjxp decay in 77^ resonance region. In full (black) 
line, we show only the resonant contribution while dotted (blue) and dashed (red) show 
the destructive and constructive interference with the non charm amplitude respectively. 



The contribution of the r]^ resonance in the x < (2.85 GeV)^ region can affect the non- 
resonant branching ratio, reducing it to 3.3 x 10"^, in the case of destructive interference, 
or increasing it to 4.2 x 10~^ in the opposite case. 

In the treatment of the decay B~ — > (j)(j)K~ , due to the complexity of the problem, 
there are uncertainties which might be important. The simplest possible approach which 
will give us a reasonable estimate of the decay rates could be the use of factorization 
model for the weak vertices and the creation of the final state by the exchange of resonant 
states. Both assumptions bring in uncertainties themselves. The model should be tested 
when more experimental data on other B decays into two vector and one pseudoscalar 
states will be available. The additional errors come from: the lack of understanding of the 
B~ — > K~ri' decay amplitude within the factorization approximation; the treatment of 
the two gluon exchange in the amplitudes of the rjK, rj'K modes and the assumptions 
on B^ 77(1490)/^" decay mechanism. The other input parameters might introduce 
about 10% uncertainty. Since the ry(1490) state gives important contribution to the rate, 
the theoretical ignorance of its coupling is potentially dangerous. 

In conclusion, we have constructed a model, based on the naive factorization and the 
exchange of intermediate resonances, with the aim to understand the decay mechanism in 
the 5^ — s> (f)(j)K~ decay with the 00 invariant mass in the region below charm threshold. 
We have found that the largest contribution in the rate comes from the decay chain B~ 
r]{rj' ,r] {1490)) K~ (j)(f)K~. Although this dominant contribution comes from the tree- 
level and penguin operators, we find that effects of the tree amplitudes are negligible. The 
interference effects of the rjc resonance with the non-resonant contribution in the region 
of the phase space with the invariant mass of the 00 state in the region (2.94 GeV)^ < 
X < (3.02 GeV)^ might decrease (or increase) the rate by ~ 10%, depending on the sign 
of the interference term. 
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